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We have attempted an approach to the chiral phase tran- 
sition of QCD using the linear sigma model as an effective 
theory. In order to get some insight into how the phase tran- 
sition could proceed, we have calculated the finite tempera- 
ture effective potential of this model in the Hartree and large 
N approximations using the Cornwall-Jackiw-Tomboulis for- 
malism of composite operators. 



I. INTRODUCTION 

It is well known the chiral symmetry of QCD is broken 
due to the small quark masses. A large amount of the 
current research activity is devoted to the chiral phase 
transition of QCD by using effective models or lattice 
techniques. An account on the physics of thermal QCD 
and further references is given in |Q . In studies of sym- 
metry restoring phase transitions the finite temperature 
effective potential is an important and popular theoreti- 
cal tool ||||. The finite temperature effective potential 
V(4>, T) is defined through an effective action T(0) which 
is the generating functional of the one particle irreducible 
graphs and it has the meaning of the free energy density 
of the system under consideration. 

A generalised version is the effective potential V(<p, G) 
for composite operators introduced by Cornwall, Jackiw 
and Tomboulis (CJT) [Q and originally was written at 
zero temperature but it has been extended to finite tem- 
perature by Amelino-Camelia and Pi || . In this case the 
effective action r(<^>, G) is the generating functional of the 
two particle irreducible vacuum graphs and the effective 
potential depends on <p(x) - a possible expectation value 
of a quantum field (f>(x) - and on G{x, y) - a possible ex- 
pectation value of the time ordered product of the field 
operator T$(x)$(y). 

Physical solutions demand that the effective potential 
should satisfy the stationarity requirmcnts 



"super daisy" graphs - figs, lb, lc respectively 
the usual tree level propagators. 



using 



dV((j>, G) 
d4> 



and 



dVjfr G) 
dG 







(1) 



The conventional effective potential results as V(4>) = 
V(4>;Go) at the solution G{<f) = Go(<fr) of the second 
equation. 

There is an advantage in using the CJT method to cal- 
culate the effective potential in Hartree approximation. 
As it was shown in a recent investigation of the \<fi the- 
ory ||, we need to evaluate only one type of graphs - 
those in fig. la - using an ansatz for a "dressed propaga- 
tor" , instead of summing the infinite class of "daisy" and 





FIG. 1. The double bubble and examples of daisy and su- 
perdaisy diagrams. 

As in many other cases in physics, in order to deal 
with the chiral phase transition of QCD we can use effec- 
tive models to describe the physical situation. The linear 
sigma model || is such a theory since it exhibits the cor- 
rect chiral properties and has attracted much attention 
recently, especially in studies of Disoriented Chiral Con- 
densates [QHIO]. ^he generalised version of the linear 
sigma model - without fermions - is called the O(N) or 
vector model and is based on a set of N real scalar fields. 
The O(N) model Lagrangian can be written as 

C = \{d^f im 2 * 2 - ^A* 4 - ea , (2) 

and in order our notation to be consistent with applica- 
tions on pion phenomenology we can identify the <&i with 
the a field and the remaining N — 1 components as the 
pion fields, so = (cr, 7Ti,... , ttjv-i)- The last term — ea 
has been introduced in order to generate the observed 
masses of the pions. 

The contact with phenomenology is obtained by con- 
sidering the case N = 4. Then the model consists of four 
scalar fields, the sigma a field and the usual three pion 
fields 7T°, tt^~ which form a four vector (cr, 7Tj) i = 1, 2, 3. 
The a field can be used to represent the quark conden- 
sate, the order parameter for the chiral phase transition, 
since they exhibit the same behaviour under chiral trans- 
formations The pions are very light particles and 
can be considered approximately as massless Goldstone 
bosons. In order to be consistent with the observed pion 
mass of w 138 McV and the usually adopted sigma 
mass of m a ~ 600 McV we choose e = j^m^, where 
,/Vr = 93 MeV is the pion decay constant. The cou- 
pling constant A and the negative mass parameter m 2 
of the model are defined as A = 3(m 2 — to 2 )// 2 and 
-m 2 = (m 2 - 3m 2 )/2 > 0. 
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II. THE HARTREE APPROXIMATION 

We examine the 0(4) model first by setting N = 4 
in the Lagrangian (Q). In the chiral limit when the pi- 
ons are considered as massless the last term is ignored. 
The four types of bubble diagrams relevant for Hartree 
approximation are given in fig. 2. 
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FIG. 2. Graphs which contribute to the effective potential 
for the 0(4) linear sigma model in the Hatree approximation. 
The numbers show the weight of each type of bubble in the 
expression for the effective potential. 



Following the procedure described in [[To| the thermal 
effective potential reads as 



V{cj>,M) = l -m 2 4> 2 + ^A0 4 + \ 



\nG-\k) 



1 



lnG- 1 (fc) + - / [V-\k)G a (k)-l] 
z Jp L J p 

+ lJ[V-\k)GAk)~3] 

2 r _ 2 

G*{k) +15 A Jc Ak) 
G a {k) f G^k)] . (3) 
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In our calculations we use the imaginary time formalism 
and for simplicity we use the following notation 



d 4 k 



d 3 k 



(2tt)4 (3^ J (2tt)3 " J0 

The tree propagators appearing above are given by 

T>- l (<P;k) = k 2 + m 2 + ^Xc/) 2 

V- 1 (^k)^k 2 + m 2 + h^ 2 , 
6 



(4) 



(5) 



and as in H we adopt the following form for the dressed 
propagator 



G-} n (4>;k) = k 2 + M 2 ^ 



(6) 



introducing a thermal effective mass M a /^. Minimizing 
the potential with respect to G a / Tr we obtain a system 
of gap equations for the thermal effective masses, which 
after performing the Matsubara frequency sums appears 
as 



Ml = m 2 + iA0 2 + ~J>(M CT ) + ^Fp(M T ) 

1 A 5A 

M 2 = m 2 + -\4> 2 + —Fp(M (T ) + —Fp{M w ) . (7) 
6 6 6 

and at the level of our approximation we keep only 
the finite temperature part of the resulting integrals, so 
Fp{M^ /<7 ) is given by 



„ t 2 r°° x 2 dx 



/ 2 exp(a 1 /2) _ l ' 



(8) 



with a = x 2 + y 2 and y — M v / a /T '. Minimizing with 
respect to the order parameter <fi we find one more equa- 
tion 



A 



= m i + _X4> J + -Fp(M a ) + -Fp(M 7T ) . (9) 

This system is solved numerically and the solution is 
given in fig. 3 
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FIG. 3. (a) The sigma and pion effective masses as func- 
tions of temperature, (b) Order parameter as a function of 
temperature. 

The finite temperature effective potential V(cj>, T) writ- 
ten in compact form is given by 

1/(0, M) = l -m 2 4> 2 + Aa^ 4 + Qp{M a ) 

+3Qp(M 7r )-^[Fp(M <7 )} 2 

-2±[Fp(M 7r )] 2 --Fp(M a )Fp(M 7r ) 

where Fp(M) is given by the equation (0) and 



Qp(M) = £ 



dxx In 



1 - exp(-a 1/2 ) 



(10) 



with a = x 2 + y 2 and y — M^u/T. The evolution of 
the potential with temperature is given in fig. 4 and we 
easily observe indication of a first order phase transition. 
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III. LARGE N APPROXIMATION 




FIG. 4. Evolution of the effective potential V(cf>, T) as a 
function of the order parameter 4> for several temperatures in 
steps of 2 MeV. The two minima appear as degenerate at 
T c ~ 182 MeV. The second minimum at (j> dissapears at 
a temperature T c2 ~ 187 MeV. 

When the chiral symmetry is broken minimizing the 
thermal effective potential with respect to dressed prop- 
agators we find the same set of gap equations (Q) for the 
effective masses as before, but minimizing with respect 
to the order parameter we obtain the following equation 



A 



A 



[m 2 + -\<j> 2 + -F p {M a ) + —Fp(M n )]<p -s = 0. (11) 

We have solved the system of equations (^) and ( |TT| ) nu- 
merically and the solution is given in fig. 5. 



When we consider N — 1 pion fields the Lagrangian 
of the model is given in equation (|J) and in the chiral 
limit we ignore the last term. Then the thermal effective 
potential looks like [flO 11 



N-l 



[D-\k)G w {k)-(N-l)} 
\{N 2 - 1 



N-l 
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X(N- 1) 
6iV 



G a (k) 



6N 



G a (k) 











(12) 



Afinimizing the potential with respect to the dressed 
propagators we find a set of nonlinear gap equations 
which in the large N approximation reduces to 



Ml = m 2 



Ml = m 2 



\<j) 2 
: \<P 2 



2A 
3 

2A 



F (M W ) 
-F {M n ) . 
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Minimizing the potential with respect to <f> we obtain 
dV(4>, M) 



(13) 



<t>M< = 



(14) 



> 





(a) 




theorem 


- — — —--JK 






in fig. 6. 
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which means that in this approximation the Goldstonc 
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FIG. 5. (a) Solution of the system of gap equations in the 
case when e / 0. At low temp the pions appear with the 
observed masses, (b) Evolution of the order parameter as a 
function of temperature. 

There is no longer any phase transition. We rather ob- 
serve a crossover phenomenon from the low temperature 
phase to the high temperature phase. 
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FIG. 6. (a) Solution of the system of gap equations in the 
large N approximation in the chiral limit. At low tempera- 
tures the pions appear as massless.(b) Evolution of the order 
parameter with temperature in the large N approximation. 
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In the high temperature phase when <fi = 0, the two 
equations become degenerate and can be used to calcu- 
late the critical temperature T c — VSfn ~ 161MeV. 

The presence of the symmetry breaking term in the 
Lagrangian (^J) results - after minimizing the potential - 
the same set of equations ( |l3| ) as before, but minimizing 
with respect to <fi we obtain the equation 



6 



2A 

—Fp(M 7T )]4> - 







(15) 



We have solved the system of equations (|l3|) and ( |l5| ) 
numerically and the solution is given in fig 7. 
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which are known to exhibit second order phase transi- 
tions However in the large N approximation the 
sigma contribution is ignored and this of course intro- 
duces errors when we calculate the critical temperature. 
In the N — 4 case which is closer to phenomenology we 
could probably obtain a better approximation in calcu- 
lating the transition temperature if we had considered 
all the interactions described by the Lagrangian of the 
model. We are planning some investigations in order to 
include these effects in the calculation of the effective 
potential. 

When we introduce a symmetry breaking term, eer, 
which generates the observed pion masses in the model, 
both in the Hartree and large N approximationds we 
found that there is no longer any phase transition. We 
rather observe a crossover phenomenon which confirms 
other results reported recently using the linear sigma 
model In this analysis there is also an indication 

of a first order phase transition in the chiral limit. First 
order phase transition is also reported in ]l4| , p"5| 
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FIG. 7. (a) Solution of the system of gap equations in the 
large N approximation in the case of broken chiral symme- 
try e / 0. At low temperatures the pions appear with the 
observed masses, (b) Evolution of the order parameter with 
temperature. 

As is the Hartree case we do not observe a phase tran- 
sition. We encounter the crossover phenomenon again; 
the difference being now the smoother behaviour of the 
order parameter. 



IV. CONCLUSIONS 

The calculation of the thermal effective potential of the 
linear sigma model in the Hartree and large N approxi- 
mations using the CJT formalism of composite operators 
proved to be very handy since we actually need to cal- 
culate only one type of diagram. In both cases we have 
solved numerically the system of the resulting gap equa- 
tions and found the evolution with temperature of the 
thermal effective masses. In the chiral limit, when the 
pion mass is ignored, the Hartree approximation predicts 
a first order phase transition, but in contrast to that, in 
the large N approximation we find a second order phase 
transition. This last observation seems to be in agree- 
ment with different approaches to the chiral phase tran- 
sition based on the argument that the linear sigma model 
belongs to the same universality class as other models 
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